F/G 14/5 


A0-Ai02 514 
UNCLA5SIFTF0 


MASSACHUSETTS INST OF TECH LEXINGTON LINCOLN LAB 

analysis of an electro-optical satellite observation.(U) 

JUN 81 6 J MAYER F19628-80-C-0002 

ETS-59 ESO-TR-BI-172 NL 


ESO-TR-81-172 


























ADA102514 


















MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


LINCOLN LABORATORY 


ANALYSIS OF AN ELECTRO-OPTICAL SATELLITE OBSERVATION 


C.J. MAYER 
Crimp 94 


PROJECT REPORT ETS-59 


4 JUNE 1981 


DTIC 



ELECTE 
AUG 6 1981 



B 


Approved for public release; distribution unlimited. 


LEXINGTON MASSACHUSETTS 













ABSTRACT 


^ A satellite has been observed through an optical telescope 
with a CCD camera, the electrical signal stored, and later 
analyzed with the aid of a digital computer. The unresolved 
satellite image is contained in 3 to 7 picture elements in fifty 
contiguous data frames. This allows a better understanding of 
the physical characteristics of the satellite image as viewed 
through the atmosphere by an electro-optical detector. 

This paper describes the data processing programs used to 
analyze the data, the analysis results, and implications for 
signal processing design. 
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I. INTRODUCTION 

A CCD camera has the potential of improving the performance 
of satellite detection systems by means of higher quant\3m 
efficiency, smaller pixel size, and more pixels on a focal plane 
than EBSICON cameras. An efficient satellite detection processor 
should be designed to exploit the signal characteristics of the 
camera it uses as an electro-optical detector. Selected portions 
of a digital satellite detection processor (IMC Signal Processor) 
and CCD camera are being developed for the GEODSS advanced tech¬ 
nology program. The output of a CCD camera, when observing a 
satellite, has been digitized, and then analyzed with the aid of 
a computer. This allowed a detailed study of the electrical signal 
that the processor would have to use prior to designing the 
processor. 

Satellite ATS-6 was observed with the 31” main GEODSS ETS 
telescope at the White Sands Missile Range. The CCD camera was 
at the prime focus and was operated in the stare mode with an 0.1 
second integration time. The ETS photometer measured the 
satellite's brightness as 10.5 Mv. The camera video was 
digitized to 12 bit resolution, sent to the MODCOMP computer 
at the ETS and stored on magnetic tape. This tape was brought 
back to Lexington and placed on the MODCOMP computer here. The 
tape was read by the MODCOMP and the data transmitted to the 
HP 9845 graphics computer and stored there. The HP computer 
was used for all the data analysis. 







The satellite observation is contained in fifty contiguous 
frames. Each frcime contains an array of 32 x 32 picture elements 
(pixels). The integration time for each frame was 0.1 second. 
Each pixel has the equivalent field-of-view of 2"7 x 2’.’7. 
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II. DATA PROCESSING 

The frames of data were read from the digital tape by the 

MODCOMP computer and transmitted, on an RS-232C interface^, to 

2 

the HP computer. Using this interface the HP computer was 
operated to look like a standard terminal to the MODCOMP. The 
major task of this interface was to write an HP terminal 
emulator program. This program has to read the data from the 
interface, store the data, and respond to the MODCOMP with 
standard terminal handshake commands. The program, "MOD-HP", 
listed in Appendix A did this at 9600 baud. 

To operate at this high a rate, each data word was read 
as a character string. The character strings were stored on 
HP tape cassettes as the data word was being read. The 1024 
character strings from each frame were separated by markers. 

Once all fifty frames of the character strings were stored on 
tape, the HP computer was disconnected from the MODCOMP, and 
was operated independently from then on. 

Two tape cassettes were required to store all the character 
strings. All the character strings were then read off the tapes 
and stored on one floppy disc using the "COPY" program listed 
in Appendix A. The character strings cannot be used in 
programs, but must be converted to the numbers they represent. 
The program "READ", Appendix A, read the character strings off 
the disc, converted the character strings to numbers and stored 

l"! Hewlett Packard, "Serial I/O Interface-98036A," 1977. 

2. Hewlett Packard, "I/O ROM Programming System 45B," 1979. 





these on tape. The "COPY" program was again used to transfer 
data from tape to disc. 

To insure that the data stored on the disc was accurate, 
a print-out of the data stored on the MODCOMP tapes was compared 
to the data stored on the HP disc. The data was identical in 
frame-to-frame and entry-to-entry comparisons. The first data 
entry in each frame of both the MODCOMP and HP data was the 
same, but was not a valid data value. To avoid problems in 
data processing the value of the first data entry, in each frame, 
was set to the mean of/ the data of each frame and is then called 
a valid data frame. 

Now the data is ready to be analyzed. The first data 
analysis program is "MAXAMP" listed in Appendix B. "MAXAMP” finds 
the maximum amplitude in each frame and prints the frame number, 
amplitude, and array position. A sample of "MAXAMP" output is 
shown in Figure 1. The frame number, e.g. T11V01:F stands for 
CCD-DIGITIZER test #11 (Til), valid data frame (V), frame number 
(01), and stored on floppy disc (F). The pixel array position 
is numbered from 1 to 1024 in each data frame. 

The next data analysis program listed in Appendix B is 
"FREDIS" which plots the frequency distribution of amplitudes 
for each frame. Figure 2 contains plots from two of the fifty 
frames. The amplitude scale adjusts for each frame depending on 
the maximum amplitude. The occurrence scale also adjusts 
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depending on the number of pixels, excluding the zero amplitude, 
for which there are from 350-400 pixels in each frame. 

The program "PEAKS" examines a frame of data and prints out 
the location and amplitude of any pixel whose amplitude is above 
a selected level. Figure 3 is a sample of "PEAKS" output. A 
program listing of "PEAKS" is in Appendix B. The frame number 
is given and the selected level (e.g. 19). The row, column, 
array position (as in "MAXAMP" output), and corresponding 
amplitude are listed. 

Program "DATAST", Appendix B, prints out the amplitude 
values of any frame by row and column. An asterislc is placed 
next to amplitudes above a selected amplitude. The output is 
displayed on the CRT or printed on the line printer. A sample of 
output from frame #11 is shown in Figure 4. The top line of 
each pair is array elements 1 to 16 and the second line is 
array elements 17 to 32. In this frame the satellite pixels 
have amplitudes of 182, 112, 32, and 40. An isolated pixel of 
noise with amplitude 32 is one array line above this. 

The "FRAME" program (Appendix B) allows a 13 x 11 pixel 
array around the satellite image of any frame to be displayed. 
After the array of pixels is displayed, the amplitudes are 
plotted first in column and then in row. Figure 5 is a print¬ 
out of the array of pixels around the satellite image from 
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frame #12. Figure 6 is a plot of the amplitude in the sixth 
column of Figure 5. The dashed line is a level at an amplitude 
of 19. Figure 7 is a plot of eunplitudes in the fifth row of 
Figure 5. The dashed line is also at an cunplitude of 19. 

Appendix C contains the pixels (5 x 6) surrounding the 
satellite image for frames 1 through 12. These were obtained 
using a modified version of program "FRAME". 
















III. ANALYSIS OF DATA 

Analysis of the output of the "MAXAMP" .program indicates 
that the highest amplitude of the satellite image moves between 
two adjacent pixels. Over the fifty frames, the highest 
amplitude switches 15 times between the two pixels, spends an 
equal amount of time in each and there is no periodicity to the 
switches. The range of highest amplitudes is 127 to 224 with 
a mean of 165. The mean of highest amplitude in one pixel is 
162, while the other pixel's highest amplitude mean is 169. 

Examining the output of "FREDIS" shows that in all frames 
most of the pixels have amplitudes from 0 to 17. The mean 
of the amplitudes from each of fifty frames ranges from 4.4 to 
4.9. This is strongly due to the large number of pixels (350 to 
400) in each frame with amplitudes of zero. The distributions 
of all fifty frames are similar. All frames have less than 20 
pixels cd}ove an eunplitude of 19. 

The output of "PEAKS" shows that in every frame the pixels 
with cimplitudes greater than 28 are either those of the satellite 
image or of isolated pixels. There are 532 pixels that have 
amplitudes greater than 19 in the fifty frames. Most of these 
pixels are either those of the satellite image or isolated 
pixels. However, 12 of the pixels in 10 of the fifty frames 
have an adjacent pixel greater than 19. In half of these cases 
the adjacent pixels are vertically related and in half horizontally. 
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The pixels with cunplitudes greater than 39 are all those of 
the satellite image. 

The output of "DATAST" shows that the noise is uncorrelated 
from pixel to pixel and from line to line. Eight-four percent 
of the pixels with an eimplitude greater than 8 have no adjacent 
pixel with amplitude greater than 8. Ninety-three percent of 
the pixels with eimplitude greater than 12 have no adjacent pixel 
with amplitude greater than 12. 

The output of "FRAME" gives information primarily about the 
pixels associated with the satellite image. Table 1 is a summary 
of the distribution of satellite image pixels by frames. The 
distribution of energy from the satellite image is very dynamic. 
The distribution changes significantly in periods of 0.1 sec to 
0.3 sec. In 61.5% of the freimes the energy distribution doesn't 
fit a circular Gaussian function. 













IV. SIGNAL PROCESSING CONSIDERATIONS 

After reviewing the analysis results, it is apparent that 
the satellite image is not stable over time periods as short 
as 0.2 sec. The instability of the image could result from three 
primary sources: satellite motion, telescope jitter, and atmospheric 
seeing. During the time of data collection, the ATS-6 satellite 
was a three-axis stabilized geosynchronous satellite. Its 
reflected light vs. time had been measured by several observers 
at the ETS and was found to be steady throughout the night, as 
would be expected from this type of satellite. Therefore satellite 
motion didn't account for the image instability observed in the 
data. The fifty frames of data were carefully scrutinized to 
detect any indication of telescope jitter. Telescope jitter 
would be detected in the data, if the pixels which contain the 
satellite image shifted together in the same direction. There is 
no observable correlation between the maximum eimplitude pixel shift 
direction and the other satellite image pixels' shift direction. 
Calculating the centroid of the satellite pixels shows that the 
centroid movement is in random directions rather than in one or 
two orthogonal directions that would characterize telescope jitter. 
Atmospheric seeing is therefore the leading candidate as the 
source of the image instability. Although few measurements of 
atmospheric seeing with time constants as short as 0.1-0.2 sec, 
and with two-dimensional spatial resolution as fine as 2'.'7 have 
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3 

been reported, one reference does report the characteristic 
image fluctuation seen in this data as atmospheric seeing. Many 
other reports on atmospheric seeing describe image fluctuation 
similar to that observed in this data. Atmospheric seeing will 
affect stellar images in the same fashion. A signal processing 
system that removes stars by comparing a star field at one time 
with the same star field at a time 0.2 sec later or longer, must 
incorporate a technique to reduce false alarms due to the star 
image instability. These false alarms are located adjacent to 
star images, and a guard band around each star that inhibits 
detections would remove these false alarms. 

Although the resolution of the telescope is 1'.'0 or smaller 
and the pixel resolution is 2'.'74, the satellite image is in 
3 to 7 pixels. Most of the energy is in 2 to 3 pixels, however 
in 74% of the frames the energy is in 4 to 5 pixels. In more 
than half the frames the distribution of satellite energy 
doesn't fit into a circular function. There is little correla¬ 
tion from frame-to-frcune of the spatial distribution of satellite 
energy. This requires that predictions of satellite SNR and 
system detection performance be modified to incorporate these 
physical realities. 


3. MIKESELL, A.H. et al, "The Scintillation of Starlight", JOSA, 
41, 689-695, 1951. 
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1 


The star and satellite images are in several adjacent pixels 
(at a 2V74 pixel size), while the noise is only one pixel long. 
This is an excellent discriminant between the satellite and high 
amplitude noise that competes with the satellite signal. If 
a threshold level of 28 was used with the pixels in the fifty 
frames, the threshold crossings would be due to the satellite 
(adjacent threshold crossings) or to the noise (isolated thres¬ 
hold crossings). The isolated threshold crossings could then be 
removed and only the satellite threshold crossings would remain. 
As the threshold is lowered, noise threshold crossings will be 
more frequent and false alarms due to noise will start to re¬ 
appear. At a threshold of 28, 100% of the noise false alarms 
are removed; at a threshold of 12, 93%; and at a threshold of 
8, 84%. The remaining false alarms due to noise can be removed 
by additional processing (streak, slope, velocity, etc.). 
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V. 


CONCLUSIONS 


The analysis of a satellite image as seen by a CCD camera 
mounted on the 31" prime focus at the ETS has been accomplished. 
The complexity of the image as a function of spatial distribution 
and time can be clearly seen. The seeing effects were worse than 
what had previously been expected at the ETS. The utility of 
having a digitized, stored signal from a sensor before designing 
a signal processor that uses that particular sensor, is obvious. 

The pixel size of a camera is important to the detection 
system; as too small a pixel will cause the satellite energy to 
be spread over many pixels and thus reduce the SNR, or if the 
pixel is too big the night sky will contribute additional noise 
and also lower the SNR. The image size depends on the atmospheric 
seeing, therefore, the selection of pixel size depends on the 
seeing expected at a site. The ideal way to measure seeing at 
the ETS would be to use an automated version of what was done 
in this report. The microprocessor based PVDS system would be 
ideal for this purpose and could provide a short exposure 
measure of seeing. Using a CCD camera mounted on the cassegrain 
position of the 31" ETS telescope would provide even more detail 
of the image intensity distribution. Binary stars could be 
viewed to eliminate some observing uncertainties and data could 
be collected on a regular basis. 
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APPENDIX A 
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0(lT=VAL<.D»a>C2,P<l.<]) 

D=1 

GOTO 290 
0<, 1 >=-l 
J1 = ..T+1 

P < J1> =POS < D»' n C P(J > +1],CHR*<13> > 

IF PCJ1>=0 THEN 480 
P<J1>=P<Jl>+P<J) 

L=LEN(:D*f: I ) CP< J> + 1, P< J1 > ] > 

IF L<4 THEN 460 
Lf =P0S';:D*< I > , CHR*< 10> > 

IF Lf<>0 THEN 390 
G=1 

GOTO 430 
G = 2 

D*a >[p<, j::'+G,P< J1 > ] = TPIM$<D#( n IP-: J>+G, P’:: J1 > ]> 
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1 

4 Id 

L =LEN<. D» < I > [ P C J > + G, P < J 1 > ] > 



4 2 m 

IF L 6 THEN 460 



4 30 

0 ( J 1 > = VflL ( D» < I ) C P >: J > ♦G, P < J 1 > ] > 



440 

ri = D+l 



450 

GOTO 470 



460 

0<; J1 '=-1 



4ro 

GOTO 500 




s=ri 



490 

GOTO 510 


1 i 

500 

NEXT J • 

NEXT J 

1 

510 

35 = 0 



520 

C5 = St +1-Si 



530 

C f =St +s 


1 

540 

ll = 0 


i 

550 

IF ';0<1 :) = -!> HND <3 = 0) THEN 670 


t 

560 

F1JR K. = Cs TO L f' ? 

PACK VARIABLES FOR OUTPUT ■: k ■ 

f 

5ro 

L =K-St ^Ss 



5S0 

IF OvL)=-l THEN 600 


1 

590 

GOTO 640 



600 

L = L+1 



6 1 0 

Cf = C f-1 



620 

Ss = Si-'-l 



6 3 0 

GOTO 580 



640 

F':k:'=o<l:- 



6 5 0 

NEXT K ! 

NEXT K 


660 

St =St +S 



670 

NEXT I ' 

NEXT 1 


6S0 

FOR M=1 TO 32 ! 

OUTPUT DATA <M> 


690 

FOR N=t1 TO 102 3+M STEP 32 



700 

IMRGE *,4D 



710 

PRINT USING 700;F';N> 



72£'I 

NEXT N 



730 

PRINT LIN(1> 



740 

NEXT t1 ! 

NEXT M 


750 

ON END #2 GOTO 790 



760 

CREATE W2»,l,,9000 



770 

ASSIGN #2 TO W2$ 



780 

PRINT #2;f(*;> 



790 

800 

"READ" (Cont.) 




10 FOR 1=10 TO 50 

! "COPY" 



20 I* = 'v'flL*':; I ) 

30 T$="T11D" 

! FRAME NUMBER 

> 


40 Plt=":T15" 

! TAPE DRIVE 



50 P2»=":F" 

60 Fl» = TtS^I»&Pl* 

70 F2»=T»&:I*&P2» 

80 PRINT F1*,F2» 

90 COPY FU TO F2* 

100 NEXT I 

110 END 

! DISC 



"COPY" 


'■ 


24 


i 

1 

i 

. 
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10 

PRINT "FRAME"," ".“POSITION" 

! "MAXAMP" i 

20 

FOR 1=10 TO 50 


30 

DIM Fa024) 

i’ 

40 

I<=VAL»<I> 


50 

VS="T11V"&I*«.":F" 


60 

ASSIGN «1 TO V$ 


70 

ON END #1 GOTO 100 


80 

READ #l;F<*> 


90 

Max=0 


100 

FOR M=1 TO 32 

! OUTPUT DATA CM> 1 

110 

FOR N=M TO 1023+M STEP 32 

1 

120 

IF F<N>>Max THEN 140 


130 

GOTO 160 


140 

Max = F<N:) 


150 

Nax=N 

j 

1 60 

NEXT N 

1 

170 

NEXT M 

NEXT M 

1 80 

PRINT V$,"MAXIMUM AMPLITUDE",Max,Nax 

j 

190 

PRINT LIN<1) 

1 

200 

NEXT I 

1 

210 

END 

1 





10 PRINTER IS 16 ! “FREDIS" 

£0 PRINT PAGE;" Thii program allows the user to enter a frame of data and ha 
"e a" 

30 PRINT " frequency distribution computed and drawn." 

40 PRINT LIN<1);" After the distribution, a beep will sound. The chart will 


50 PRINT "remain on the screen and the program will halt until you press CONT 

It 

60 PRINT "This will give you time to digest the information on the screen." 

70 PRINT LIN'^O;" If you want the output of this program dumped from the CRT 

II 


30 PRINT "to the internal thermal printer, then answer Y to the question below 


90 PRINT "If you don'’t want hardcopy answer N. Then press CONT." 

100 INPUT "Do you want hardcopy ‘:Y/'N>?",At 

110 IF UPCS'.:A$> = "Y" then Yhard 

120 IF UPC*<A»>="N" THEN Nhard 

130 BEEP 

140 GOTO 100 

150 Yhard: Hardcopy=i 

160 GOTO 130 

170 Nhard: Hardcopy=0 

130 N=224 

190 CALL Dr i wer <N, Har dcopy :• 

200 END 

210 SUB Driuer<N,Hardcopy> 

220 RANDOMIZE 

230 DIM A(:N>,L$<:N>, Uni ts$C30],V< 1024 > 

"FREDIS" 


25 
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I'40 INPUT "Which fruit,e 01 to 50?“,Pr$ 

250 V^ = "TllV"&Fr•*^^":F" 

2>;0 ASSIGN #l TO V$ 

270 REAIi 
2S0 Ma:=V';i> 

2'ji0 FOP J=1 TO 1024 
500 Ma;:=MAX(Max, VC J> * 

510 NEXT J 
520 MAT A=ZEP 
550 High=0 
540 Lhigh=409fc" 

350 FOP K=1 TO 1024 

5 6 0 IF C y C K '> > High.' AND <■ V C h > . L h i g h > THEN 3 S 0 
370 GOTO 390 

3 S £' H 1 g h = V i, K > 

390 NEXT K 

400 FOR Jl=l TO 1024 

410 IF VC Jl ::.=High THEN 4 30 

420 GOTO 440 

430 A C H i gh ) =A <. H i gh • + 1 

440 NEXT J1 

450 IF High=0 THEN 490 

4fc'0 Lhigh=High 

470 High=0 

4S0 GOTO 350 

490 PRINTER IS 0 

500 PRINT 

510 PRINT V$ 

520 PRINT "MAX AMP" 

550 PRINT Max 
540 PRINTER IS 16 

550 FOR 1=20 TO Max STEP 20 ! ENTER THE TITLES 

560 L* C I :> =VAL* C I > 

570 NEXT I 

5S0 Unit£*="Number of Pixeii" 

590 Ti11 €#="Amp 1itudes" 

600 N = M ax 

610 CALL Bar C AC*> , L* c * :;i, N, Urn t s*, Ti 1 1 €*;> 

620 BEEP 

6 5 0 PAUSE 

640 IF Hard':opy=l THEN DUMP GRAPHICS 
650 EXIT GRAPHICS 
6 6 0 S U E E XI T 

670 SUE Bar (A(■ * >, LJ c * >, N, Uni t i$, Ti 1 1 e$> 

680 PLOTTER IS "GRAPHICS" 

690 GCLEAP 
700 GRAPHICS 
710 PEN 1 
720 Surii=0 


"FREDIS" (Cont.) 
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7 50 l'1ax = fl' 1 ' 

740 FOR 1=1 TO N 
750 Sijni=Sc4rii + fl( I > 

760 Max = MRX'.'Max, |4'; I > > 

770 NEXT I 
780 Lv=LGTtMax> 

780 SETGU 

300 LOCRTE 15,120,40,95 • Set screen boundaries and scale 

810 CLIP 0,123,40,100 
820 MOVE 61.5,97 
830 LORG 5 

840 LRBEL USING 850;Tit1e« 

850 IMAGE #,K 

360 SCALE 0, N+. 5,0, Max+Max 10 

870 Xtic=0 ! No tic marks on horizontal axis 

380 Testyt i c=FRACT (Ly > + <;Ly<0) 

390 Vt i c = 10 ■" < INT CLy > - 1 > *< 1 + 1.5*< <Testyt i c >. 39794) AND < Test yt i c 69397 + 4+ C <; Te 

St ytic >=.69897 ) AND < Test ytic < = .87506))+6.5*< Test ytic >.87506)> 

900 ! Compute the uertical tic marks 

910 CALL LaxesCXtic,Ytic,0,0,1,1,2,0,N+.5,0,Max+Max/10) 

920 LORG 4 
930 FIXED 1 
940 FOR 1=0 TO N 
950 MOVE I-.5,0 

960 IPLOT 0,A<I),-1 ! Draw the bar for the ith data item 

970 IPLOT 1,0 
980 IPLOT 0,-A<I) 

990 MOVE I,A<I) 

1000 NEXT I 
1010 DEG 
1020 IMAGE K 

1030 LDIR 90 ! Rotate the label 90 degrees 

1040 SETGU 
1050 LORG 5 
1060 MOVE 5,70 

1070 LABEL USING 1020;Units$ ! Write the display units being 

1080 ! used along the uertical axis 

1090 CLIP 0,123,0,100 

1100 SETUU 

1110 FOR 1=1 TO N 


"FREDIS" (Cont.) 
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f 


1120 MOVE I,-Max 20 
1130 LORG 8 
1140 PEN 1 

1150 LDIR 90 ! Rotate the labels 90 degree 

1180 LABEL USING 1020;L*<I> 

1170 NEXT I 
1180 BEEP 
1190 SUBEXIT 

1200 SUB Laxes < Xt 1 c , Yt i c , Xorg, Yorg, Xmaj , Ymaj , M j rit i c s J re, Xoi i n, Xmax, Ym i r, Vmax 
1210 IF <: Xm i ri > = Xmax > OR Ytn i n > = Ymax ) THEN SUBEXIT 
1220 Xfudge=.02*<Xmax-Xmin> 

1230 Yfudge=.02* < Ymax-Ymin> 

1240 AXES Xtic,Ytic,Xorg,Yorg,Xmaj,Ymaj,Mt nticsire 
1250 DEG 

1260 IF NOT Xtic THEN Lately 

1270 LabeIxlIF SGN(Xtic)»-l THEN Parx 

1280 LDIR 90 

1290 LORG 8 

1300 GOTO 1330 

1310 Parx; LDIR 0 

1320 LORG 6 

1330 FOR I=Xorg TO Xmax STEP ABS<Xtic) 

1340 MOVE I,Yorg-Yfodge 
1350 LABEL USING 1360;I 
1360 IMAGE «,K 
1370 NEXT I 

1380 Lately: IF NOT Ytic THEN SUBEXIT 

1390 IF SGN<Ytic>*-l THEN Pary 

1400 LDIR 0 

1410 LORG 8 

1420 GOTO 1450 

1430 Pary: LDIR -90 

1440 LORG 6 

1450 FOR I = Yorg TO Ymax STEP ABS<Ytic;> 

1460 MOVE Xorg-Xfudge,I 
1470 LABEL USING 1360;I 
1480 NEXT I 
1490 SUBEXIT 


"FREDIS" (Cont.) 
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10 

DIM F<:i024> 

! "PEAKS" 

20 

F*="T11V22:F" 


30 

L=19 


40 

PRINTER IS 0 


50 

PRINT F», "PERKS ABOVE",L 


60 

PRINT 


70 

PRINT "ROW","COLUMN","POSITION", 

"AMPLITUDE" 

80 

PRINTER IS 16 


90 

ASSIGN #1 TO F$ 


100 

ON END #1 GOTO 120 


1 10 

READ #1;F<*> 


120 

FOR M=1 TO 32 

! OUTPUT DATA 

1 30 

C = 0 


140 

FOR N=M TO 1023+M STEP 32 


150 

CsC + 1 


160 

IF F<N)>L THEN 200 


170 

IMAGE «,4D," " 


180 

PRINT USING 170;F<N> 


190 

GOTO 250 


200 

IMAGE #,4D,“*" 


210 

PRINT USING 200;F<N) 


220 

PRINTER IS 0 


230 

PRINT M,C,N,F<N) 


240 

PRINTER IS 16 


250 

NEXT N 


260 

PRINT LIN<1) 


270 

NEXT M 

! NEXT M 

280 

PRINTER IS 0 


290 

PRINT LINO) 


300 

END 



"PEAKS” 
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4 


10 DIM Fa024> ! "DflTHST" 

20 INPUT "AMPLITUDES HIGHER THEN",L 

30 FOR J=10 TO 50 

40 PRINT LINO) 

50 PRINT "FRAME NUMBER",J 

60 PRINT LINO) 

70 S«=VAL*0) 

80 ASSIGN #l TO "T1 1V“«<S$8,":F" 

80 ON END #1 GOTO 110 

100 READ #l;F<*) 

110 FOR M=1 TO 32 ! OUTPUT DATA (M) 

120 FOR N=M TO 1023+M STEP 32 

130 IF FCN)>L THEN 170 

140 IMAGE #,4D," " 

150 PRINT USING 140,'FCN) 

160 GOTO 180 

170 IMAGE «,4D, 

180 PRINT USING 170;F<N) 

190 NEXT N 

200 PRINT LIN(l) 

210 NEXT M ! NEXT M 

220 NEXT J 

230 END 


•DATA5T" 


30 
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1C1 INPUT “LEVEL^L 

20 FOR 1=10 TO 50 

30 DIM F< 1024:j 

40 I» = VFlL»<:i> 

50 V*="T11V"&I»4":F" 

60 0SSIGH il TO Vt 

?0 ON END #1 GOTO 90 

80 READ #l;F';*> 

90 GRAPHICS 

100 LINE TYPE 1 

110 SCALE 0,12,0,15 

120 FOR M=15 TO 25 

130 C=0 

140 FOR N=M TO 1023+M STEP 32 

150 C=C+1 

160 IF FCNXL THEN 200 
170 MOVE M-15,C 

180 LABEL F<N) 

190 IF C=13 THEN 210 

200 NEXT N 

210 NEXT M 

220 WAIT 3000 

230 GCLEAR 

240 SCALE 0,13,0,230 

250 FOR M=19 TO 22 ! PLOT COLUMNS 

260 WAIT 200 

270 C=0 

280 GCLEAR 

290 MOVE 0,L ! PLOT LEVEL 

300 LINE TYPE 4 

310 DRAW 13,L 


! “FRAME" 


! DATA FROM DISC 


! LABEL AMPLITUDE 
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320 LINE TYPE 1 

330 MOVE 1,150 

340 LfIBEL "COLUMN" 

350 move 0,0 

360 FOR N=M TO 1023+M STEP 32 

370 C=C+1 

380 PLOT C-1,F-;N> 

390 PLOT C,F<N) 

400 NEXT N 

410 WAIT 500 

420 move 0,0 

430 NEXT M 

440 GCLEAR I plot ROWS 

450 move 0,0 

460 FOR Nl=239 TO 335 STEP 32 ^ LEVEL 

470 MOVE 0,L 

480 LINE TYPE E 

4’^0 DRAW 13,L 

500 LINE TYPE 1 

510 MOVE 1,150 

520 LABEL "ROW" 

530 MOVE 0,0 
540 FOR M=1 TO 13 
550 N=N1+M 
560 PLOT M-1,F<H> 

570 PLOT M,F<:N> 

580 NEXT M 

590 move 0,0 

600 WAIT 500 

610 GCLEAR 
620 NEXT N1 

630 NEXT I 

640 END 
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172 
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